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Abstract 

Using a waveguide-type sideband-separating receiver (2SB receiver) on the Tokyo- 
NRO 60-cm telescope (renamed the AMANOGAWA telescope), we carried out simul- 
taneous observations in the 12 CO( J = 2 — 1) and 13 CO( J = 2 — 1) lines over the Galactic 
plane / = 10° - 245° along 6 = 0° with a 3^75 grid. Using the 12 CO( J = 1-0) data 
of Dame et al. (2001), who used a beam size almost the same as ours, we show 
12 CO(J = 2 - l)/ 12 CO(J =1-0) and 13 CO(J = 2 - l)/ 12 CO(J = 2-1) intensity 
ratios on the I — v map and the intensity correlations among the 12 CO(J = 2 — 1), 
13 CO( J = 2 — 1), and 12 CO( J =1 — 0) lines. As a result, a linear correlation between 
12 CO(J=1-0) and 12 CO(J = 2-l) and a curve correlation between 12 CO(J = 2-l) 
and 13 CO (J = 2 — 1), as produced by most of the data, are found. We investigate 
these correlations with simple radiative transfer equations to ascertain a number of 
restrictions on the physical quantities of molecular gas on a galactic scale. 
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1. Introduction 



In order to understand the dynamics and evolution of galaxies, it is important to study 
the distribution and the physical condition of interstellar gases. In particular, molecular gas 
is important because it forms stars, which occupy about 90% of the visible mass of their host 
galaxies. The Schmidt law gives the relation between the star formation rate and the mass 
surface density of the gas (Schmidt, 1958). As a result of observational investigation, its power 
index has been determined to be about 1.4 (Kenicutt, 1998). 

In one empirical scenario for the formation of massive stars, stars are formed in a molec- 
ular core whose density is remarkably high, at about 10 7 cm -3 (e.g. Gracia-Carpio et al. 2008; 
Muraoka et al., 2009). On the other hand, the fact that the Schmidt law remains valid for a 
giant molecular cloud (GMC) scale of ~ 10 2 pc suggests a relationship between the molecular 
cores and the gas properties averaged over a GMC. In order to investigate this relationship, it 
is necessary to obtain high resolution data on a galactic scale. The Milky Way Galaxy (MWG) 
is the best target source for this purpose, because it can be observed with the highest spatial 
resolution and S/N ratio. We can elucidate the relationship between the physical condition of 
GMCs and the formation of stars in their environment through detailed observation. If the data 
cover the greater part of the MWG and does not omit the less dense gases, we can compare 
them with the data of external galaxies, whose data are smoothed over a large area with a 
telescope's beam. 

Dame et al. (2001) carried out a survey of the entire MWG in 12 CO(J =1 — 0) and 
showed the distribution of molecular gases, including less dense gases. Observations coupled 
with this data make it possible for us to investigate the physical conditions of the molecular 
gases. If the observations are carried out in the line of another transition of CO or its isotope, 
we can minimize the uncertainty introduced by chemical abundance. 

Sanders et al. (1986) and Jackson et al. (2006) carried out surveys that covered a large 
area of the MWG in 12 CO( J =1 — 0) and 13 CO( J =1—0), respectively. However, the coverage 
of their observations was limited to a part of the first Galactic quadrant and they did not 
observe the higher- J transitions of the same molecular species. 

Sakamoto et al. (1995) carried out a survey in 12 CO( J = 2 — 1), whose coverage was also 
limited to / = 20° — 60°, using the same beam size as Cohen et al. (1986), and they found a 
decrease in the 12 CO( J = 2-1) / 12 CO( J =1-0) intensity ratio with an increase in the Galacto- 
centric distance (Sakamoto et al., 1997). They could not, however, offer a detailed discussion 
of the 2-dimensional distribution due to the insufficient sensitivity of their receiver. 

We upgraded the Tokyo-NRO 60-cm telescope used by Sakamoto et al. (1995) by in- 
stalling a waveguide-type sideband-separating (2SB) SIS receiver. It then had a sensitivity 
that was 10 times better, and the two frequency bands could be observed simultaneously 
(Nakajima et al., 2007). Using this system, we began a new survey of both 12 CO(J = 2 — 1) 
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and 13 C0( J = 2 — 1) over the Galactic plane. 

As the first result of this survey, in this paper we present the data on the Galactic 
equator {b = 0°). In addition, using the survey data, we discuss the condition of the "typical" 
molecular gas in the MWG without a complicated molecular excitation model, as the first step 
in our investigation. 

In Section 2, we describe the specifics of our observations. The data and comparisons 
of the three CO lines are given in Section 3. Using these results, in Section 4, we discuss 
the physical condition of the typical molecular gas in the Galactic disk. Finally, we present a 
summary in Section 5. 

2. Observations 

2.1. Instrumentation 

The telescope we used is the Tokyo-NRO 60-cm telescope (VST-1), which was specifically 
designed for rapid sky surveys on the 200 GHz band. It was built in 1990 and upgraded in 2006 
by replacement of the receiver system. In this paper, we refer to the new telescope system as 
the AMANOGAWA telescope, which is an acronym for 'A Miniature Antenna at NObeyama 
for GAlactic objects and the milky WAy'. 

A 2SB receiver developed by Nakajima et al. (2007) was installed on our telescope. This 
is the first operated receiver of this type used for scientific observations on the 200 GHz band. It 
enables us to obtain two frequency band datasets on both sides of the local oscillator frequency; 
the upper sideband (USB) and the lower sideband (LSB). In addition, its waveguide-based 
design produces a lower noise temperature. The typical system noise temperature, including the 
atmosphere, is about 200 K at the zenith. Our spectra were obtained using two 2048-ch Acousto- 
Optical Spectrometers (AOS) developed by Kaiden et al. (2006). Their frequency coverage 
and resolution are 250 MHz and 250 kHz, respectively. These values correspond to a velocity 
coverage of 326 km s _1 and a velocity resolution of 0.3 km s _1 at 230 GHz, respectively. The 
basic characteristics and detailed specifications of the AMANOGAWA telescope are reported 
in Nakajima et al. (2007). 

Although Nakajima et al. (2007) used the new moon, we found that it was not suitable 
for measurements of the main beam efficiency and the half power beam width (HPBW) of our 
telescope. The measured antenna temperature of the new moon changed from month to month, 
although the intensities of our calibration source were stable (see subsection 2.3). This can be 
attributed to the sensitive change of the beam averaged brightness by the moon phase and 
the position on the moon. The measured antenna temperature also was not stable for the full 
moon. Therefore, we measured the efficiency and HPBW using the sun in January 2009. In 
order to avoid saturation of the SIS mixer, we installed a wire grid in front of the receiver. 

The sun is larger than the main beam of the telescope and the sidelobe level is sufficiently 
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low. Measurement of the sidelobe level showed it to be less than -18 dB relative to the peak. 
In this case, the main beam efficiency ^ m b can be estimated as 

_ /mainbcam P^Q ^ _ /sun -Pd£l _ ^A,sun /-, \ 

^ mb " k,P^ ~ VSUU ~ h«P^ ~ T sun ' [l) 
where P is the beam pattern, T£ sun is the measured antenna temperature of the sun, and 
T sun is the intrinsic brightness temperature of the sun. The difference between rj^ and 77 sun is 
estimated to be <1 %. The T^ sun is measured using attenuation rate g of the wire grid as 

J A,sun — rp J-atmy 

gj-hot — gJ-sky 

where T no t is the temperature of the hot load and T s k y is the brightness temperature of blank 
sky. The temperature of the atmosphere T atm is assumed to be equal to T hot . The estimated ?7 m b 
is free from the value of g. We carried out the measurement several times using the different 
values of g that ranged from 0.25 to 0.75. Using a T sun of 6256 K (USB) and 6299 K (LSB) 
(Linsky, 1973), we obtained a r/ mb value of 0.725±0.029 for the USB and 0.767±0.038 for the 
LSB. 

The HPBW was measured at 8f7±0f4, a value that is quite close to that of the 1.2 m 
telescope used in the CfA CO( J =1 — 0) survey conducted by Dame et al. (2001). The spatial 
resolution corresponds to 20 pc at 8 kpc from the sun. 

The dataset we used in this work was obtained in 2007 and the measurement of the beam 
was revised in 2009. We maintained the same telescope optics, however, during this period. 
This means that the revised ?7 m b and HPBW are valid for all the data in this paper. 

2.2. Observations and Data Reduction 

The observations were carried out in the first season of the AMANOGAWA-2SB Galactic 
Plane Survey from January to April, 2007. We observed the 12 CO(J = 2-l) and 13 CO(J = 2-l) 
lines simultaneously. The sky coverage was / = 10° — 245° along 6 = 0°, which is the entire 
range of longitude accessible from the observation site, except near the Galactic center. The 
grid spacing was 3 .'75 or 0.4 times that of HPBW. The doppler tracking velocity was fixed at 
ulsr = km s _1 for each line. We obtained the data of the velocity range from -160 km s^ 1 to 
+160 km s -1 in all spectra. 

All data were taken in a position-switching mode, in which off positions about 5° from 
the Galactic equator were selected from those given by Dame et al. (2001). After a 10-second 
integration period, the on position was switched to the off position. We summed up 6 to 8 sets 
to obtain the final integration data. 

Due to the stable receiver system and offset optics, the baselines of the spectra are quite 
flat, so that only linear baselines were subtracted. After this procedure, Gaussian smoothing 
was applied to all the profile data in order to reduce the noise level. For the 12 CO data, 
a Gaussian with 1.3 km s~ x velocity width was applied, which is the same as the velocity 
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resolution of Dame et al. (2001). The resultant noise level was a tms ~ 0.05 K. For the 13 CO 
data, a Gaussian with a width of 2.6 km s _1 was applied to improve the noise level for the 
weaker emission features. The resultant noise level was <r rms ~ 0.04 K. 

2. 3. Calibrations 

To confirm the accuracy and internal consistency of these data, we monitored the ref- 
erence points in Ml 7, W3, and OriKL every 2-3 hours in an observation day. The measured 
intensities of these sources are consistent within ±7% and the velocity accuracy is better than 
0.3 km s" 1 for each span in which the receiver tuning was maintained (Yoda, 2008). However, 
the apparent intensities of the calibration sources changed slightly following each receiver tun- 
ing. This apparent change is likely due to differences in the image rejection ratio (IRR) of the 
receiver. In an ideal case, the IRR would be infinity and the receiver line intensity would be 
at a maximum. Therefore, we adjusted the intensity scale in each tuning span so that it would 
match with the maximum that had ever been measured. Since the IRR of both sidebands 
changes independently, we determined the scaling factor individually. After this correction, the 
main beam temperature at a reference point in OriKL, (/, b) = (209.°028765, — 19.°375472) was 
35.51±3.44 K in 12 CO(J = 2-1) and 10.76±0.84 K in 13 CO(J = 2-1). These values were 
derived using 5-point fitting around the peak position (Yoda, 2008) to reduce the pointing error 
effect. 

The optical and radio pointing observations (see Nakajima et al., 2007) were made every 
1-2 months, and we revised the pointing parameters so that the pointing accuracy would be 
better than 1'. We checked this using the data of the reference points shown above, and 
confirmed that it would be about V. 

3. Results 

3.1 . Overall View of Data 

Although the observed area has one pixel extent along the Galactic latitude, it had many 
pixels along the LSR velocity. Therefore, an overall view of our data can be shown in longitude- 
velocity (/ — v ) maps. Figure 1 and Figure 2 are I — v maps of the main beam temperature 
in 12 CO(J = 2 — 1), T 12 , and in 13 CO(J = 2 — 1), T 13 , respectively. The velocity widths of a 
pixel on the I — v plane are 1.3 km s _1 in the 12 CO line and 2.6 km s _1 in the 13 CO line. For 
comparison, we also show an I — v map of the main beam temperature in 12 CO( J =1 — 0), Ti-o, 
which was constructed based on Dame et al. (2001) (Fig. 3). The velocity width is 1.3 km s~ x 
and the grid spacing is 7.5 arcmin, which is twice as sparse as ours. 

Since the 12 CO( J = 2—1) and 13 CO( J = 2 — 1) datasets were obtained simultaneously 
using the same receiver horn of the same telescope, there are no relative pointing errors, and 
the difference in beam size is less than 5%. Moreover, our sampling and beam size coincide 
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with those in the observations of Dame et al. (2001). We are therefore able to compare these 
data directly without any image processing. 

The l — v maps of T12, T13 and Ti_o look very similar. They are not, however, completely 
proportional to one another. For example, the T12 intensity in the region of / = 10° — 20° and 
ulsr > 60 km s _1 is much weaker than that of Ti-o, compared with its intensity in other regions. 
We put the 12 CO(J = 2 - l)/ 12 CO(J = 1-0) line intensity ratio, #21/10, on the I — v plane in 
Figure 4 and the 13 CO(J = 2 — l)/ 12 CO(J = 2 — 1) ratio, #13/12, on the same plane in Figure 5. 
The clipping levels are > 5cr rms for both lines. Both figures are drawn with a velocity width of 
5.2 km s _1 to reduce the noise level. This resolution is narrow enough to permit investigation 
of the Galactic structure because the typical velocity width of a cloud is similar to or wider 
than the resolution. 

Figure 4 shows that most of the data pixels have a similar value. However, the standard 
deviation of #21/10, 0.147, is larger than was expected from the system noise. #21/10 actually 
changes from place to place and it shows a difference in the physical condition of the molecular 
gas. The average value of #21/10 is 0.613. In limited regions around the terminal velocity, or 4 
and 6 kpc from the Galactic center, for example (I, wlsr) — (50°, 60 km s _1 ) and (25°, 110 km 
s" 1 ), some pixels show a remarkably high #21/10 of about 1.0. 

Figure 5 shows that the distribution of #13/12 roughly corresponds to that of the line 
intensities of each line. For example, the peaks of T 13 at (/, i>lsr) — (31°, 100 km s" 1 ) and (54°, 
20 km s _1 ) correspond to the peaks of #13/12- There seems to be a general trend in which the 
more intense pixel shows the higher #13/12. The average value of #13/12 is 0.166. 

3.2. Correlations between Line Intensities 

In order to discover the characteristics of the Galactic disk gas, the correlations between 
the intensities in several lines are useful. We therefore constructed an intensity correlation plot 
between T12 and T13 in Figure 6 and between Ti_o and T12 in Figure 7. In both figures, we only 
use the data in I = 10° — 65° in order to reduce the effects of local molecular clouds and the 
Cygnus region. This approach is valid for investigating the global trend of the MWG, because 
this longitudinal range accounts for 84% of the mass of the molecular gas in the entire observed 
area. 

In the correlation between T12 and T13, most data points are aligned along a curve with 
increasing #13/12 (Fig. 6). This means that #13/12 is higher for the more intense data point, as 
we mentioned in the previous subsection. We will discuss the interpretation of this correlation 
curve in the next section. 

In the correlation between T12 and Ti_o, most data points are aligned in a linear corre- 
lation (Fig. 7). This means that most of the gas shows only a small variation in #21/10- The 
slope of the regression line is 0.640. The observational error for the slope due to statistical 
noise, baseline error, and intensity fluctuation is estimated to be ±0.051. The dispersion of all 
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Fig. 1. The I — v map of the 12 CO(J = 2—1) line intensity. Each panel shows a different longitudinal 
range shown as the abscissa labels. The colors are shown in log-scale at the point where T12 > 3<7 rms . 
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Fig. 2. The Z — w map of the 13 CO( J = 2 — 1) line intensity. The data for T13 are displayed in the same 
way as in Figure 1, though the color coding is different. 
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Fig. 3. The I — v map of the 12 C0( J = 1 — 0) line intensity produced from the observation data of Dame 
et al. (2001). The panels are displayed in the same way as in Figure 1 and 2, though the color coding is 
different. 
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Fig. 4. The I — v distribution of the ratio of T12 to Ti_o- Both lines are averaged in a 5.2 km s _1 velocity 
bin. The color coding is shown at the top-left. The data with T12 < 5tr rms ~ 0.15 K or T1-0 < 5er rms ~ 0.2 
K are shown in black. 
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Fig. 5. The data of X13 and X12 are averaged in a 5.2 km s _1 velocity bin. The ratio between them, 
^i3/i2i is displayed on the I — v plane. The data with T13 < 5er rms ~ 0.15 K or T±2 < 5cr rms ~ 0.15 K are 
shown in black. 



11 



2 












12 3 

Intensity in '-CO I / M 



4 



Fig. 6. The intensity correlation between T13 and T12. Each datum is averaged in a 5.2 km s 1 velocity 
bin. The data plot forms a curve with increasing slope. 

the measured data points gives ±0.058 for the slope, which includes the intrinsic dispersion of 



In Figure 7, a number of data points are located above the regression line. This suggests 
that there may be a trend between line intensity and #21/10 ■ Next, we constructed a plot on 
the T12 — #21/10 plane for the data on the I — v map with T12 ^ 0.25 K (Figure 8). Note that 
the error of the ordinate is not uniform along the abscissa, as shown in the lower panel of Fig. 
8. The error of the ratio is smaller than ±0.05 for T 12 ^1 K, while it reaches ±0.1 for T 12 ^ 0.5 
K. The error of T12 is almost uniform at a typical value of <r rms ~ 0.05 K. 

At T 12 ^1 K, the data points are scattered from #21/10 —0.4 to 1.2. The scattered #21/10 
is real, because it is much larger than the observational error. In contrast, for T12 ^3 K, #21/10 
converges to about 0.6. 

4. Discussion 

We found that there are correlations between the line intensities in 12 CO (J = 2 — 1) and 
13 CO(J = 2-1) and those in 12 CO( J = 2-1) and 12 CO( J = 1-0), although some data show 
deviations beyond the observational error. These relationships should reflect some constraint 
on the physical condition of the molecular gas in the Galactic disk. In this section, we discuss 
only the physical condition of the greater part of the gas that displays these relationships, or 
the global trend. We will not offer detailed discussion of the gas with a large deviation, which 
we will address in forthcoming papers. 

To ascertain the actual relationship between the parameters of the physical condition 
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Fig. 7. The intensity correlation between Ti_o and Ti 2 . Each datum is averaged in a 5.2 km s 1 velocity 
bin. The data plots forms a straight line. 
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Fig. 8. The upper panel (a) shows the correlation between Tyi and i?2i/io- Each datum is averaged in a 
5.2 km s _1 velocity bin. The lower panel (b) shows the error of i?2i/i0i Ai? 21 / 10 . The Ai? 21 / 10 depends on 

the value of T12, Ti_o and their noise levels through Ai? 2 i/io =^/ jr^ 2 AT 1 2 2 + ^ J - 2 AT 1 2 _ . The average 

value of Ai? 21 / 10 in a T 12 bin is represented by a fillcd-in circle and its dispersion due to Ti_ 0l AT 12 and 
ATi_o are represented by bars. 
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of the molecular gas, we adopted a means of deriving them without the use of gas excitation 
models such as the large velocity gradient (LVG) model (Goldreich & Kwan, 1974). We will 
give relationships empirically from the results presented in the previous section, based on the 
simple radiative transfer equation with the beam dilution effect: 

T = r,T c (l -exp(-r)), (3) 

where 77 is the beam filling factor and r is the optical depth in the line. Since our observation 
beam is larger than the typical size of a molecular cloud, the beam filling factor strongly affects 
the line intensity. T c is the Rayleigh- Jeans line temperature 1 that is introduced to simplify the 
expression. 

In subsection 4.1 and 4.2, we derive the empirical relationships between the parameters 
of the molecular gas. Based on these, we discuss the typical physical condition of the gas 
in subsection 4.3. In subsection 4.4, we briefly discuss the data points that deviated from 
the typical relationship presented in the previous section. In this paper, we assume that the 
emitting region of a beam is uniform. 

4-1. Correlation between T 12 and T 13 and Two Empirical Relationships of the CO gas 

From the correlation curve between the line intensities in 12 CO (J = 2 — 1) and 13 CO (J = 
2 — 1), we can obtain two restrictions on these lines. One is that the optical depth in the 12 CO 
line is much larger than that in the 13 CO line. If it is not, then the correlation curve must be 
close to a straight line. This is consistent with the relative abundance of these two molecular 
species; it is reported to be about 60 in the Galactic disk (e.g., about 59 by Liszt & Lucas, 
1998; about 76 by Stahl et al., 2008). With this abundance, the optical depth in the 12 CO line 
is expected to be several tens of times as thick as that in the 13 CO line. 

This means that the 12 CO lines should be thick in most cases. Although many observa- 
tions have already revealed that the 12 CO line is optically thick in the nearby molecular clouds, 
our findings show that it is also thick in most of the molecular gas beyond the boundaries of 
mapped molecular clouds. 

In this case, we can obtain the following relationships, using equation 3: 

Tiz = 7713^,13(1 -exp(-ri 3 )), (5) 
T12 = 7712^12 • (6) 




The Rayleigh- Jeans line temperature T c is defined as 

(4) 

where v and T ox are the frequency and the excitation temperature of the transition, respectively. Tb g is 
the brightness temperature of the cosmic microwave background radiation (CMB) in the Rayleigh- Jeans 
approximation. Since the CMB is the 2.725 K blackbody, T hgA2 = 0.19 K, T bg ,i3 = °- 22 K and T b g ,i-o = 0.84 
K at the corresponding frequencies. T c is close to T ox in the case that T ox 3> ^ and T ox 3> Xt, g . 
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Fig. 9. The distribution of x 2 given by equation 9 on the a — f3 plane, which are the two parameters of 
the model expressed in equations 7 and 8. The \ 2 value is low along a line and takes its minimum at (a,/3) 
= (0.4, 1.3K). 

The approximation in the second equation is valid within 5%, if Tu > 3. 

We should note that some data should be optically thin in 12 CO(J = 2 — 1). In particular, 
the data points with #21/10 ^1-0 are difficult to explain if the line is opticaly thick. Such data 
points, however, are less than 7% and such gas is beyond the scope of the present discussion. 

Equations 5 and 6 have five parameters, but the effective number of free parameters 
in only three, because both r] 13 T c 13 and 7712^,12 affect the equations as single parameters. 
Therefore, the correlation curve between the 12 CO and the 13 CO lines gives a number of rela- 
tionships among 77i 3 T C)13 , rj 12 T Cl i2 and n 3 . 

Within the first order approximation, they must be 



where a and (3 are constants. To estimate these two cnstants, we evaluate \ 2 as defined by 



where the superscript "obs" represents the observed value, and "exp" represents the value 
expected from the observed value using equation 7 and 8. The calculation is carried out for the 
averaged T° 3 bs over each T° 2 bs bin with a 0.1 K step in order to avoid over- weighting the data 
points with small T° 2 bs and to obtain a globally fitted curve. 

The resulting distribution of \ 2 on f ne a ~~ ft plane is shown in Figure 9. It indicates 
that the x 2 value is small along a line. Although (a, /3) = (0.4, 1.3 K) gives the smallest x 2 , 
the correlation curve is also produced in the case of (0.3, 0.7 K) or (0.5, 2.8 K) (see Fig. 10). 



— = — = a, 

?7l2-7c,12 

^13^c,13 = P ■ T13, 



(7) 
(8) 




(9) 



What is the astrophysical meaning of these two constants? Equations 7 and 8 would 
represent the "common sense" of the molecular gas in the Galactic disk. In the following 
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Fig. 10. The correlation curve produced by our model equations 7 and 8. The blue, cyan and green solid 
lines are correspond to (a, (3) = (0.4, 1.3 K), (0.3, 0.7 K), and (0.5, 2.8 K), respectively. The red dots are 
the T12 — T±s intensity correlation used in the \ 2 fitting. The fillcd-in circles in each color correspond to 
Ti 3 =0.01, 0.30, 0.50, and 1.00, respectively. 

discussion, we elaborate this idea using the terminology of astrophysics. 

The constant a is the ratio of 7713^,13 an d Vi2T c , 12- We divide it as a product of 7713/7712 
and T c>13 /T c>12 . 

7713/7712 is the ratio of the projected sizes of the emitting region. As shown in the first 
part of this subsection, the 12 CO line is optically much thicker than the 13 CO line. This suggests 
that the emitting region in the 12 CO line should be larger than that in the 13 CO line. Sawada 
et al. (2001) estimated that in the Galactic center region, 7713/7712 is about 2/3 or 0.7. Although 
the actual value of the 7713/7712 in the Galactic disk should be different from that in the Galactic 
center region, it is difficult to estimate. For now, then, we therefore assume that it is about 
0.7. 

In the case that 7713/7712 is constant over the Galactic disk, the constant a means that 
T Ci i 3 /T c i2 is constant. This does not directly indicate that T CX;1 3/T eX) i 2 , because T c is not 
linear to T cx . For 7713/7712 = 0.7, a = 0.4 ±0.1 gives T c l3 /T c i2 = 0.6 ±0.1, which corresponds to 
T ex ,i 3 /T CXi i2~0.7±0.1. 

Since the 12 CO line can be emitted from a less dense gas than the 13 CO line can, 7712 
should be larger than 7713 for a cloud with an envelope of less dense gas, and 7712 should be equal 
to 7713 for a cloud without the envelope. Therefore, the upper limit of 7713/7712 should be unity. 
In this case, T eXi i 3 /T eXi i2 is about 0.5. The lower limit of 7713/7712 is difficult to estimate without 
a model of the density structure of a cloud. If we presume 7/13/7/12 = 0.5, then T eXi i 3 /T eXi i2 is 
about 0.8. 
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The other constant should show the condition of the 13 CO emission, because it is 
the linear coefficient of 713 to 7713^13. The derived value of f3 changes from 0.7 K to 2.8 K 
depending on the value of a. Equation 8, which defines /?, is qualitatively reasonable once we 
consider the following discussion. A larger gas cloud gives a larger r and a larger 77. In a gas 
cloud of fixed size in which the gas is subthermally excited and in which photon trapping is 
essential for molecule excitation, a denser gas gives a larger r, and a larger r gives a higher 
T cx and T c for the cloud. Although it is difficult to show which is more affected by T13, 7713 or 
T c ,i3, equation 8 may indicate a relationship between T13 and 7713, because both parameters are 
directly affected by cloud size. 

4-2. Ratio of T^and Ti_oand Excitation of 12 CO 

From the linear correlation between the line intensities in 12 CO (J = 2 — 1) and 12 CO (J = 
1 — 0), we can obtain a restriction on 12 CO excitation. We showed #21/10 — 0.640 ±0.058 in 
Section 3.1. Since the 12 CO lines are optically thick, #21/10 is expressed as 

#21/10 - ~ • (10) 

r ?l-0 J c,l-0 

Although we cannot evaluate the exact value of the beam filling factors, 7712 should be close to 
771-0, because the excitation conditions for both emission lines are similar and the apparent size 
in both lines are close for a model cloud (Sakamoto et al., 1997). In this case, we can deduce 

#21/10 = J"' 12 ■ (11) 

J c,l-0 

Therefore, #21/10 = 0.640 ± 0.058 indicates that 12 CO( J = 2 — 1) is subthermally excited. 

In the previous subsection, we concluded that 12 CO( J = 2 — 1) is optically thick. This 
large optical depth allows for the small probability that the emission line photons may escape 
and make the gas thermalized. However, moderate opacity can give both subthermal excitation 
and optically thick approximation for the line intensities. For example, the gas with Tk ^ 20 K, 
n(H 2 )= 10 2 - 0_2 - 5 cm~ 3 and iV( 12 CO)/d?j = 10 16-17 cm -2 km -1 s shows subthermality and 7~i 2 — 5, 
using the LVG model. 

A constant T C)12 /T Cjl „o does not always give a constant T ex . 12 /T CXjl _o, because T c is 
linear to T cx only when T cx ^> hu/k and T ex 3> Tb g . Figure 12 shows the loci of T Ci i2/T Ci i_o = 
0.640 ± 0.058 on the T cx l2 — T CXi i 2 /T CX) i_ plane. The value of T eXtl2 /T eXjl _ increases near 
T ex ,i2 = T^g. From our data, we found that most of the molecular gas shows almost constant 
#21/10- If T cx ,i2 is not fairly uniform over the Galactic disk, Fig. 12 suggests that there is 
a lower limit on it, because the gas that is colder than this temperature must have a higher 
T , eX) i2/T , C x,i-o than the warmer gas in order to keep the same #21/10 value. There is a non- 
negligible amount of molecular gas that is warmer than 20 K. Therefore, we can conclude that 
T eXt i2 is warmer than 13 K. In this case, T eX) i2/T ex> i_o is 0.7 and T eXj i_o is higher than about 19 
K. 
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Fig. 11. Regression line derived by the \ 2 fitting, which has a slope of 0.640, and lines with a slope of 
0.640±0.058 (dashed lines). The background is the Ti_o — T1.2 intensity correlation plot shown in Figure 
7. 

Combined with T CXjl3 /T eXil 2 given in the previous subsection, we can estimate the lower 
limit of 13. For T eXil3 /T exl2 ~ 0.7, it is about 9 K. Although we did not make observations in 
13 CO( J =1 — 0) with the same beam size, its excitation temperature should be higher than 12 
K, if the ratio of the excitation temperature of 13 CO( J = 2 — 1) and 13 CO( J =1 — 0) is similar 
to T eX) i 2 /T eXj i_o, — 0.75. Using higher resolution (44 arcsec) data for the molecular clouds 
in the 4kpc-ring region, Rathborne et al. (2009) estimated that the excitation temperature 
in 13 CO(J =1 — 0) line is 9 K. Our estimation is not inconsistent with theirs, because these 
estimations may have large ambiguity due to the many assumptions involved. 

4-3. Gas Kinematic Temperature and Typical Gas Density. 

As shown in the previous section, the lower limit of T ex l _ is about 19 K. This means 
that the lower limit of the gas kinematic temperature Tk should also be about 19 K. We believe 
that Tk is not much higher than T ex> i_o, since 12 CO( J = 1 — 0) is optically thick and the energy 
level of J = 1 is only 5.5 K greater than the ground level. 

It is very difficult to estimate the gas density in molecular clouds from low resolution 
data that is insufficient to resolve individual clouds. We can, however, discuss the typical gas 
density in a cloud with the help of the derived values given in previous subsections. 

With low resolution data, the beam filling factor is of essential importance in deriving 
the physical property of interstellar gas. In the optically thick case, it was given as the ratio 
of the measured brightness temperature to the Rayleigh- Jeans line temperature, as shown in 
equation 6. Since all of data points show T12 < 7 K and about 80% of the data points in our 
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Fig. 12. The theoretical loci of fl 2 i/io = 777^ = 0.640 (solid line), 0.640-0.058 (dahsed line), and 
0.640+0.058 (dotted line) on the T cx .i2 — t^'^o pl &ne J as calculated with equation 4. If T cx ,i2 changes 
between 13 K and 40 K, the constant T ox ,i2/r o x,i-o gives the same -R21/10 within the observational error. 
If T eX[ i2 is colder than 13 K, the constant -R21/10 requires a higher T eXi i2/T eX) i_o for the colder gas. 

sample show T 12 < 1.5 K, 7/ 12 must be less than 0.7 and should be less than 0.14 for most of 
the data points, though some of them may be different from those for gas with a homogeneous 
property, which is discussed in this paper. 

We can estimate the volume filling factor ( which is the volume fraction of line emitting 
region to the sampling volume in an observation beam, from the beam filling factor rj, when 
the distribution of the emitting region is isotropic in the sampling volume. We can neglect the 
overlapping of clouds in the sky because the actual beam filling factor is much smaller than 
unity. Thus we obtain 

C~^i, (12) 

where r is the spatial resolution at the subcentral point, or at the average distance to the two 
positions corresponding to the same velocity, given by the beam size, and s is the path length 
of a data cell. Using equation 6, we estimated the upper limit of 7712 for each data point. Using 
a kinematic model of the MWG, we can also estimate r and s for each data point. 

Many investigations of the Galactic plane and external galaxies use the assumption that 
Xi_o is proportional to the H 2 column density iV(H2) through the CO-to-H 2 conversion factor 
X for the first order approximation. Using the path length s, we can estimate the H 2 volume 
density averaged over the beam and the velocity bin as 

(n(H 2 )) beam = — Ti_ . (13) 

s 

According to the definition of the volume filling factor, (n(H 2 ))b ea m is equal to the gas density 
in molecular clouds (^(H 2 )) c i ou d diluted by £■ Therefore, we can estimate the lower limit of 
(n(H 2 )) doud for each data point. A histogram of this shown in Figure 13 with X = 1.8 x 10 20 
cm -2 (K km s" 1 )" 1 (Dame et al., 2001). The estimated lower limit of (n(H 2 )) c i ou d is typically 
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Fig. 13. Histogram of the lower limit of the gas density in molecular clouds (n(H2)) c ioud as calculated by 
(n(H2))bcam/C The histogram shows a peak at (n(H2)} c i ud —100cm -3 . 

more than 100 cm" 3 , which is consistent with the results derived from the data of a nearby 
molecular cloud (e.g., Sakamoto et al., 1994) or molecular gas in a Galactic disk of a galaxy 
using a CO excitation model (e.g., Muraoka et al., 2007). 

4-4- Data Points that Deviate from Main Trend 

Until the previous subsection, our discussion focused on the global trend of the data. A 
number of data points in Figure 6-8, however, deviate from the trend. In this subsection, we 
make observations on the characteristics of these data points. 

In subsection 3.2, we showed that -R21/10 has a large intrinsic dispersion at a weaker T 12 
(Fig. 8). We believe that such data points with weak 7\ 2 are located in the inter-arm region, 
because the spiral arm is traced by concentrations of molecular gas, which indicates intense 
T 12 . In fact, observations of external galaxies have shown that intense features in a CO line 
are distributed on the spiral arm (e.g., Kuno et al., 2007). The data points with intense T12 
are continuous in the l — v plane and they correspond to spiral arm (Fig. 1). It is certain 
that #21/10 is controlled by the temperature and/or density, though it is difficult to factorize 
their respective contributions. Therefore, the large dispersion at T 12 suggests that the physical 
condition is more heterogeneous in the inter-arm region than in the spiral arm. 

In Figure 7, a number of plots deviated from the regression line toward the higher T12. 
These data points are also located far from the majority in the T 12 — -R21/10 plane (Fig. 8). One 
would not expect the data points in an area of Ti 2 ~ 1 - 4 K and R 21/1Q > 0.86 to be optically 
thick, because #21/10 ~ 0.86 requires T exl2 > T cxl „ in the optically thick case. On the / — v 
plane, most of them are located near Hll regions either at the terminal velocity in / ~ 25° — 30° 
and I ~ 50° or near W49 and W44. 
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The T12 — T13 correlation plot (Fig. 6) also has a number of data points that are located 
far from the correlation curve. Some plots at Tyi - 5 K showed remarkably lower T13 than the 
global trend. These data points show a fairly higher -R21/10 than the averaged value and should 
have lower opacity in 12 CO( J = 2 — 1). They correspond to a molecular cloud at I ~ 13.5° and 
vlsr ~ 10 km s -1 , which is overlaid on a low mass X-ray binary, GX13+1 (White et al., 1978). 
This cloud would exist under the unique physical condition. High resolution observations would 
be helpful to resolve this issue. 

5. Summary 

In order to investigate the global physical condition of the MWG, we carried out simul- 
taneous observations of 12 CO and 13 CO J = 2 — 1. We covered the Galactic plane at b = 0° over 
I = 10° - 245° with a 3.75 grid using the AMANOGAWA telescope. Using the 12 CO( J =1-0) 
dataset of Dame et al. (2001), which has almost the same spatial resolution, we showed the 
distribution of the 12 CO( J = 2- l)/ 12 CO( J = 1-0) and the 13 CO( J = 2- l)/ 12 CO ( J = 2 - 1) 
ratios on the I — v plane and the intensity correlations between their intensities. 

We found that for the majority of data points, the intensity correlation between 
12 CO( J = 2 - 1) and 13 CO( J = 2 - 1) is aligned along a curve, and that between 12 CO( J = 2 - 1) 
and 12 CO(J = 1 — 0) lies along a straight line with an intensity ratio of -R21/10 = 0.640. From 
these results, we conclude that the optical depth in the 12 CO (J = 2 — 1) line is much greater 
than that in 13 CO( J = 2—1), and most of the data points follow the two simple relationships 
shown in equations 7 and 8. Using the best fit values of parameters a and /3 in these equations, 
we estimated typical values for the physical condition of the molecular gas in the Galactic disk 
as follows: the excitation temperatures in 12 CO( J = 1 - 0) , 12 CO( J = 2 - 1) , and 13 CO( J = 2 - 1) 
are higher than 19 K, 13 K, and 9 K, respectively. The gas kinetic temperature is higher than 
19 K. The H_2 volume density in a cloud is more than 100 cm -3 . The beam filling factor is 
much less than 0.7. 

We found a number of data points that deviate from the global trend. This means that 
the gas temperature and density in these locations are different from the typical values we 
obtained. Some of these are located close to Hll regions or an X-ray source, while the others 
are located in the interarm region. 

The AMANOGAWA telescope is a joint project of the University of Tokyo, the 
National Astronomical Observatory of Japan, Osaka Prefecture University, and Tokyo Gakugei 
University. We especially thank the staff at the Nobeyama Radio Observatory, Nobeyama Solar 
Radio Observatory, and the Advanced Technology Center for their support in operating the 
telescope. We would also like to thank the builders of the original telescope system, before the 
upgrade. Finally, we thank our referee for the valuable advice given to improve this paper. This 
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